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ABSTRACT 

We present a series of smoothed particle hydrodynamical models of GO. 253+0. 016 (also known as 
"The Brick"), a very dense molecular cloud that lies close to the Galactic Centre. We explore how 
its gas and dust temperatures react as we vary the strength of both the interstellar radiation field 
(ISRF) and the cosmic ray ionisation rate (CRIR). The cloud has an extent in the plane of the sky of 
roughly 3.4 pc x 9.4 pc. As its size along the line-of-sight is unknown, we consider two cases. In our 
fiducial, high-density model, we adopt a depth along the line-of-sight of 3.4 pc, and in the low-density 
model, we assume an extent along the line-of-sight of 17 pc. To recover the observed gas and dust 
temperatures, we find that the ISRF must be around 1000 times the solar neighbourhood value, and 
the CRIR must be roughly 10 _14 s _1 , regardless of the geometries studied. For such high values of 
the CRIR, we find that cooling in the cloud's interior is dominated by neutral oxygen, in contrast to 
standard molecular clouds, which at the same densities are mainly cooled via CO. Our results suggest 
that the conditions near GO. 253+0. 016 are more extreme than those generally accepted for the inner 
500 pc of the galaxy. 
Subject headings: stars: formation 



1. INTRODUCTION 

The environmental conditions in the Galactic Cen- 
tre (GC) provide an extreme test of our c urrent under- 
standing of the star formation p r ocess (e.g. | Pap adopou- 
Tosl[20T0l IKrumholz et al.||2012| |Longmore et al.| |2013; 
Krui jssen et al.||2013| ). With both stronger background 
radiation fields and higher cosmic-ray fluxes compared 
to clouds in the solar neighborhood, star formation is 
predicted to occur at higher volume and column densi- 
ties than is typical in a standard giant molecular cloud 
dElmegreen et al.|2008[ ). 

One notable example is GO. 253+0. 016 (also referred to 
as the "The Brick"), which displays both extremely high 
column an d volume densities, yet very little sign of star 
formation (Glisten et al.|1981 ; Lis et all 1994 [Longmore 



et al. 2012). Despite the current lack of star formation, 
the physical conditions in this object are thought to be 
similar to thos e required for the form ation of massive 
stellar clusters ( |Longmore et al.||2012 ). 

In this paper we investigate the influence of the ex- 
treme GC environment on the thermodynamics of dense 
and massive molecular clouds, in an attempt to better 
understand the initial conditions for star formation in 
the inner molecular zone. We adopt values for the inter- 
stellar radiation field and the cosmic ray ionisation rate 
that are significantly higher than those measured in solar 
neighborhood molecular clouds. For more fundamental 
parameters such as the mass, dimensions, and turbulent 
velocity dispersion of the c louds, we take th e value s for 
GO. 253+0.016 reported by |Longmore et al | ( [2012] ). In 
contrast to the other clouds in the GC, the apparent 
lack of star formation in GO. 253+0. 016 makes it an ideal 
candidate for studying the effects of the environmental 
conditions on the thermal balance of the cloud. 




Fig. 1. — Column density, and mean gas and dust temperatures 
in our fiducial cloud setup (simulation '1' in Table fTJ, with the 
ISRF set at 1000 G , and the CRIR at 3 x 1Q- I4 s- 1 . 



2. COMPUTATIONAL METHOD 
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TABLE 1 

Summary of the simulations 



Model 


2 T 

[pc] 


3 T 

[pc] 


4 U 
[pc] 


^mirijO 

[cm" 2 ] 


6 n 
[cm -3 ] 


7 IlSRF 

[Go] 




8 Icr 
Is" 1 ] 




9 x(R 2 ) 


10 x(CO) 




n x(C+) 




12 x(0) 




1 


9.4 


3.4 


3.4 


3.6 x 10 23 


3.5 x 10 4 


1000 


3 


x 10" 


14 


0.477 


1.5 x 10- 


5 


1.33 x 10" 


5 


3.06 x 10" 


-4 


2 


9.4 


3.4 


3.4 


3.6 x 10 23 


3.5 x 10 4 


100 


3 


x 10" 


15 


0.500 


9.59 x 10" 


5 


8.02 x 10" 


6 


2.24 x 10" 


-4 


3 


9.4 


3.4 


3.4 


3.6 x 10 23 


3.5 x 10 4 


1000 


3 


x 10" 


16 


0.500 


1.11 x 10" 


4 


3.08 x 10" 


6 


2.01 x 10" 


-4 


4 


9.4 


3.4 


17.0 


7.3 x 10 22 


6.7 x 10 3 


100 


3 


x 10- 


16 


0.500 


6.65 x 10" 


5 


1.76 x 10" 


5 


2.53 x 10" 


-4 


5 


9.4 


3.4 


17.0 


7.3 x 10 22 


6.7 x 10 3 


100 


3 


x 10- 


15 


0.496 


2.62 x 10- 


5 


2.41 x 10- 


5 


2.94 x 10" 


-4 


6 


9.4 


3.4 


17.0 


7.3 x 10 22 


6.7 x 10 3 


1000 


3 


x 10- 


16 


0.497 


6.10 x 10- 


5 


2.52 x 10- 


5 


2.59 x 10" 


-4 



Note. — 2,3,4 Initial physical dimensions of the cloud. 5 Minimum column density, measured along the shortest axis. 6 Initial 
hydrogen nuclei number density. 7 Strength of the interstellar radiation field, in units of the local value. 8 Cosmic ray ionisation 
rate. sU * 11 *^ p ma i fractional chemical abundances in the cloud, measured at the point at which the first core goes into runaway 
collapse. These are quoted with respect to the number of H nuclei. A fully molecular gas therefore has x(H2) = 0.5. The total 
carbon and oxygen abundances in the models are 1.4 x 10 — 4 and 3.2 x 10 — 4 respectively. 



We perform our simulations using the smoothe d par 



tide hydrodynamics (SPH) code Gadget2 (Springel 
2005). We have modified the code to include time 
dependent chemistry and a treatment of the main heat- 
ing and cooling processes (described below). We have 
also in cluded an implementation of t he TreeCol algo- 
rithm (Clark, Glover fc Klessen||2012 ), to obtain column 



density maps ol the sky as seen by each SPH particle. 
These maps (including total, H 2 and CO column densi- 
ties) are used to calculate the influence of the interstellar 
radiation field (ISRF) on the gas and the dust. 

We assume for simplici t y tha t the spectral sha pe of 
the ISRF follows " 



Draine 



^ fl978l ) in the UV and |Black| 
( 1994 ) at longer wavelengths. We denote the solar neigh- 
bourhood value of the strength of the ISRF as Go, 
and perform simulations with field strengths 100 Go and 
1000 Go (see Table [lj. Note that this multiplicative scal- 
ing is done equally at all wavelengths. For our dust 
model, we use a combination of values from Ossenkopf & 
Henning|(|1994|) ( non-coagulated, thick ice mantle gr ains) 
lor w avelengt hs longer th an 1 //m, and from Mat his, 
Mezger fc Panagia] ([1983 ) at shorter wavelengths. To 
compute the visual extinction, we use the relationship 
Am = 5.348 x 10- 22 (7V H ,tot/l cm" 2 ), where 7Vh 



V = o.o^o a ±u V iV H,tot 
the total hydrogen column density 



Drake 1978 Draine & Bertoldi 199oT For simplicity, 



jlBohlin, Savage 



we do not account tor any changes in the extinction 
curve that may occur due to dust coagulation. For the 
cosmic-ray ionisation rate (CRIR), we adopt a value of 
I CR n = 3 x 10 -17 s 



1 as our solar ne ighbourhood value 
and assume that 



van der Tak fc van Dishoeck| 2000| ) 
each ionisation event deposits 20 ev of energy into the 
gas ( jGoldsmith fc L anger || 19 78] ) . The depende nce of the 
CRIR on column density is highly uncertain ( jPadovani^] 
Galli & Glassgold 2009), and we assume for simplicity 
that no attenuation occurs. We do not include the effects 
of ionization by hard X-rays, as this does not appear to be 
a major heat source in the Gal actic Center, given the rel- 
ative l y low X-ray luminosity fR odrigue z^Fernandez et al. 
2004| (Schleicher, Spaans fc Klessenj|2010| 
For the chemistry we a dopt the reduced CO net work of 



Nelson fc Langer (1999). Details can be found in Glover 
fc Clark (2012b]), and a description of how the chemistry 



interact s with the ISRF via the TreeCol algorithm is 
given in|Glover fc Clark (2012a). 



3. INITIAL CONDITIONS AND MODEL PARAMETERS 



For the initial conditions in this study, w e take the 
cloud properties derived in Longmore et al. (2012) for 
GO. 253+0.0 16 as a guide: a size ol 9.4 pc by 3.4 pc, and 
a mass of 1.3 x 1O 5 M . The clouds are simulated us- 
ing 2 x 10 T SPH particles, and so ou r mass resolution is 
M res = 0.65 M flHubber et al.|2006| ). We adopt a simple 
rectangular cuboid geometry, matching the longer of the 
two observed dimensions with the x-axis in the simula- 
tions, and the shorter with the ?/-axis, such that all the 
clouds have particles placed initially from to 9.4 pc in 
x (L x ) and to 3.4 pc in y (L y ). In the z direction we 
adopt two values for the extent of the cloud, since the 
true dimension of GO. 253+0. 016 along the line-of-sight 
is unknown. Our first choice is to make the z-axis the 
same length as the y- axis, yielding a mean hydrogen nu- 
clei number density no = 3.5 x 10 4 cm -3 . This is the 
setup used in our 'fiducial' clouds. Our second choice is 
to make z the longest axis, with L z = 17.0 pc. These 
clouds have an initial density of 6.9 x 10 3 cm -3 and are 
our 'low-density' clouds. All the clouds are given non- 
thermal support in the form of a turbulent velocity field, 
which has a power spectrum P(k) ex k~ 4 . The turbu- 
lence is permitted to decay as the cloud evolves. We 
fix the initial 3D turbu lent velocity dispersion based on 
the observational data: Longmore et al. (2012) report a 
linewidth of 15.1 kms -1 for GO. 253+0. 016, equivalent to 
a ID velocity dispersion of 6.4 kms -1 , and hence to a 3D 
velocity dispersion of 11.12 kms -1 , assuming isotropic 
turbulence. 

We perform three simulations for each of our two cloud 
models, varying the strength of the ISRF and the mag- 
nitude of the CRIR. An overview of the simulations can 
be found in Table [I] A central assumption here is that 
the shape of the radiation field and the cosmic ray en- 
ergy spectrum are the same locally and in the Galactic 
Centre, and that it is only the normalization of each that 
changes. 

In view of the high densities probed by our initial con- 
ditions, we assume that the hydrogen in our clouds starts 
as H2. However, we start with the carbon in the form 
of C + , and allow it to self-consistently evolve to form C 
and CO. As we discuss in ^5j the clouds are already in 
chemical equilibrium at the point at which we perform 
our analysis. 

4. GAS AND DUST TEMPERATURE 
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Fig. 2. — Gas (blue) and dust (red) temperatures as a function of x. The top row contains the clouds that have the fiducial setup (x is the 
longest axis), while the bottom row contains the low density clouds (those with z as the longest axis). The lines denote the mass-averaged 
temperature along the line of sight. Vertical bars denote the l-cr dispersion. 

excitation ammonia lines suggest that GO. 253+0. 016 has 
a complex gas temperature structure, with components 
up to 400 K, that has yet to be modeled (E. Mills, 2013, 
private communication). What environmental condi- 
tions are required to produce such temperatures? 

The typical features of our cloud are illustrated in 
Fig. [I] which shows the column densities of one of 
the clouds in the x-y plane (i.e. integrating along z), 
and the accompanying mean gas and dust temperature 
maps. This cloud is our most extreme case studied, with 
Iisrf = 1000 Go, and Icr = 1000Icr,o, and our 'fidu- 
cial' cloud geometry. However, the features of this cloud 
are mirrored in our other simulations - the clouds have 
a hot skin and a relatively cool interior, and are highly 
structured by the supersonic turbulence. The images in 
Fig. [I] are taken just as the first collapsing core exceeds 
a density of around 10 8 cm -3 , and so represent the state 
of the cloud at the onset of star formation. All the other 
clouds in this study will be presented at the same point 
in their evolution. 

In Figure [2j we show the gas and dust temperatures 
in the clouds as a function of the position along the x- 
axis. The most obvious feature of these profiles is that 
the gas and dust have different temperatures throughout 
the cloud. They are not thermodynamically coupled on 
the scales shown here, consistent with the observations 
mentioned above. 

The profiles also reveal how the environment affects 
the cloud temperature. We see that the cosmic rays 
are responsible for heating the gas, while the ISRF is 
primarily responsible for heating the dust. Such a re- 
sult is expected. The high column density of this cloud 
means that photo-electric emission in the cloud interior 
is strongly suppressed, as the UV photons responsible 




Fig. 3. — Gas (blue) and dust (red) temperatures as a function 
of density in our fiducial cloud (model '1' in Table 



Using Herschel observations, LongmoreetaLl (2012) 
show that the dust temperature varies smoothly from 
19 K in the cloud centre to 27 K at the edge. Ob- 
servational constraints on the gas te mperature of 
GO. 25 3+0. 016 have existed for some time: Glis ten et al.| 
(1981) derive rotation temperatures of ~45K using am- 
monia transitions, corres ponding to an average kineti c 



temperatures of ~80K (|Walmsley fc Ungerechts 



A recent formaldehyde survey ( |Ao et al.|[7013| finds 



1983). 



average kinetic temperatures of 65-70 K, which agrees 
within the uncertainties. However, observations of high- 
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for it are readily absorbed near the surface of the cloud. 
As such, the ISRF can play only a minor role in directly 
heating the gas. On the the other hand, as the cosmic 
rays have no attenuation in our model, they are free to 
heat the cloud's gaseous interior throughout. The ISRF 
can, however, heat the dust at the centre of the cloud, 
as this heating comes primarily from longer wavelength 
photons, which are able to penetrate much further than 
the UV photons. In summary, for clouds with such an 
extreme column density as GO. 253+0. 016, the heating of 
the dust and gas is effectively split into two components. 

Our 3D modelling results suggest that for our fiducial 
cloud model, the environmental parameters that best re- 
produce the observed temperatures are Iisrf = 1000 Go, 
and Icr= 1000 Icr,o- Reducing either of these values by 
a factor of ten results in gas or dust temperatures that 
are too low to agree with the observations. 

One potential source of error is simply that we have 
underestimated the extent of GO. 253+0. 016 along the 
observed line-of-sight, and so the true effective column 
of the cloud is much smaller than we are assuming in the 
fiducial models. However we find that similar environ- 
mental conditions are also required when we consider our 
lower-density version of GO. 253+0. 016. These models are 
shown on the bottom row of Fig [2] Even in these lower 
column density clouds, we see that the ISRF is mainly 
responsible for determining the dust temperatures (i.e. 
there is very little gas-dust thermodynamic coupling), 
and the CRIR is mainly responsible for determining the 
gas temperatures. Our dust temperatures are now a lit- 
tle higher than the observed values throughout the cloud, 
suggesting that for this geometry the Iisrf would need to 
be lower than 1000 Go- However we see that by 100 Go, 
the ISRF is already too low to explain the observed tem- 
peratures. Also, we see that Icr= 100 Icr,o results in 
a gas temperature of around 30 K in the interior of the 
cloud - again, this is inconsistent with the observations. 

Figure [2] also shows that the geometry of the cloud af- 
fects the temperature gradients along the cloud. This is 
particularly evident when one looks at the gas tempera- 
ture, especially when Iisrf is high (see e.g. the bottom 
right panel). This implies that it should be possible to 
constrain both the total ISRF and the cloud's geome- 
try by fitting the gradient of the gas temperature in the 
cloud modelling. Such a study is worth revisiting once 
maps of the gas temperature with sub-parsec resolution 
become available. 

Finally, we note that both the gas and dust tempera- 
tures can vary considerably along a line of sight from the 
averages shown in Fig. |2j This can already be seen in 
the images in Fig. [T] However, we also show in Fig. [3] 
how the temperatures vary as a function of density in our 
fiducial case. We see that at high densities (> 10 6 cm -3 ), 
once the dust and gas thermally couple, the temperatures 
can be relatively cold. 

5. HEATING AND COOLING PROCESSES 

In this section we investigate the heating and cool- 
ing processes for the gas in more detail. The dominant 
processes that govern the gas temperature are shown 
as functions of density in Fig. [4] for the two most ex- 
treme cases: our fiducial cloud (no = 3.5 x 10 4 cm -3 ) 
with Iisrf = 1000 Go and Icr = 1000Icr,o? and one of 
the lower-density clouds (no = 6.7 x 10 3 cm -3 ), with 




1 3 1 4 10 5 1 6 1 7 

n [cm -3 ] 



photo elec CO 

H 2 photo dis — — C + H 2 

H 2 form T DG a A GD a 

cosmic ray T pdV □ A pdV □ 

shocks 

Fig. 4. — Processes responsible for heating and cooling the gas 
in two very different cloud models (clouds 1 and 4 from Table 
Heating processes are shown in red and orange and cooling 
processes are represented in blue. Two processes - pdV work and 
gas-dust thermal coupling - can produce either heating or cooling 
depending on the circumstances. Heating and cooling associated 
with compression and expansion are denoted by T p &y and A p dy, 
respectively, while the transfer of energy from the gas to the dust 
is denoted by Aqd and that from the dust to the gas by Tdg- The 
plotted quantities represent the median values at each density. 

Iisrf = 100 G and I C r = 10I C r,o- 

In both clouds, the dominant heating processes follow 
a broadly similar pattern. At the lowest densities, which 
represent the outskirts of the clouds in these simulations, 
the dominant heat source is photoelectric emission from 
dust grains. This falls off sharply as we move to higher 
densities as a result of the increasing extinction as one 
moves into the cloud's interior. At slightly higher densi- 
ties, the heating caused by cosmic rays starts to dominate 
the thermal balance. In the case of the hotter, denser 
cloud, this process remains the main heating source until 
we reach a number density n = 10 8 cm -3 , correspond- 
ing to our resolution limit. In the lower density cloud, 
embedded in the less extreme environment, shock heat- 
ing becomes the main source of heat input to the gas at 
densities above n ~ 10 5 cm -3 . Note that since neither 
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compression nor shock-heating are dominant in the high 
CRIR case, the temperature of the cloud cannot be used 
to determine its age. 

When we compare the main cooling processes, we also 
find some similarities. In the low-density outskirts, where 
the gas is warm and there is little CO, we find that C + 
and neutral oxygen emission are the main coolants, as in 
the low-density ISM. Given the high densities and tem- 
peratures of the cloud's skin, and the fact that we start 
with the hydrogen in molecular form, we also find that 
H2 can be an effective coolant at the outskirts. 

As we move into the cloud, however, the gas tem- 
perature drops and the C + recombines to form C and 
then CO. The identity of the dominant coolant therefore 
changes. In the low-density cloud, CO cooling dominates 
in this slightly denser regime, just as is the case in local 
molecular clouds. In the denser cloud model, however, 
CO never dominates; instead, atomic oxygen becomes 
the main coolant. This difference in behaviour is a result 
of the CRIR in these two clouds. In the higher density 
cloud, the much higher CRIR creates many He + ions 
that react destructively with the CO molecules: 



CO + He+ -> C + + O + He. 



(1) 



It also keeps the gas warm enough to excite the fine struc- 
ture lines of atomic oxygen. In the lower density cloud 
with the much lower CRIR, both of these effects are less 
important, and hence atomic oxygen never becomes the 
dominant coolant. Since we need a large CRIR to ex- 
plain the observed gas temperatures, the implication is 
that the cooling of gas in GO. 253+0. 016 (and probably 
also in other Galactic Centre clouds) is dominated over 
a significant range in densities by emission from atomic 
oxygen. 

At very high densities, dust becomes the most effec- 
tive source of cooling. However this does not occur un- 
til the gas density is more than an order of magnitude 
higher than the mean cloud density, and hence we ex- 
pect that T gas = Idust only in the densest gas within 
GO. 253+0. 016, with most of the volume of the cloud hav- 
ing T gas 7^ ^dust- As already noted, this expectation is 
supported by the available observational data on the gas 
and dust temperatures. 

The effect of the clouds' environment on the chemical 
balance is summarised in Table [I] We see that strong 
ISRFs and CRIRs have little effect on the H2 fraction, 
and so we would expect the true molecular state of the 
cloud to be relatively independent of the environment. 
However, the CO fraction varies by around an order of 
magnitude in the models, implying that its ability to 
trace the molecular state of the gas is a strong function 
of the environment. Since the clouds initially have all 
of their carbon in the form of C + , one might argue that 



we have simply ended our simulations too soon to pick 
up all of the CO. However, we see that in the clouds 
with smaller CRIRs over half of the carbon is in CO, 
suggesting that there is sufficient time available for it to 
form in large quantities. As such, the low CO abundances 
in the clouds with high CRIR are due to real differences 
in their chemical evolution. 

6. DISCUSSION 

Our results suggest that the CRIR and ISRF around 
GO. 253+0. 016 should be 1000 times the solar neighbour- 
hood values, in order to obtain temperatures consistent 
with the values derived from observations. Such radia- 
tion and CR fields could be produced by enhanced star 
formation activit y, higher stellar densities, or some com- 
bination of both. [Yusuf-Zadeh et al.| ( |2009| ) measured the 
star formation rate (SFR) in the GO to be 50-100 times 
the local SFR. If the CRIR and ISRF are set solely by 
star formation, our results suggests that the local SFR 
near GO. 253+0. 016 is about an order of magnitude higher 
than the mean SFR of the central molecula r zone ( Morris 
fc Serabyn|[l996l |Yusuf-ZadehltaTp555|> . 



Similarly, the CRIR that we require is significantly 
higher than the values found for local dense clouds. How- 
ever, there is considerable observational evidence that 
the ionizatio n rate is h i gher i n the Galactic Centre. 
For example, Oka et al. (2005) estimate a value of 2- 
7xl0 ~ 15 s~ 1 in diffuse gas along s evera l GC sightlines 



while Yusuf-Zadeh, Wardle & Roy (2007) infer a value of 
2-50x 10 _i4 s _i within GC clouds, based on observations 
of the fluorescent 6.4 keV Ka iron line. Our required 
value of a few times 10 -14 s _1 is compatible with these 
values, given the large uncertainties. 

Our models also suggest that the neutral oxygen emis- 
sion coming from GO. 253+0. 016 should be significantly 
higher than that seen typical molecular clouds. This 
could provide an independent test of the models pre- 
sented in this paper. 
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